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Abstract 

We study the /3-ensemble that represents conformal blocks of Liouville theory on the sphere. This quantity is 
related through AGT conjecture to the Nekrasov instanton partition function of 4d Af = 2 SU (2) gauge theory 
with four flavors. We focus on the large N limit, equivalent to the Nekrasov-Shatashvili limit where one of 
the ^-background deformation parameters is vanishing. A quantized Seiberg-Witten differential form is defined 
perturbatively in ft as the singular part of the /3-ensemble resolvent. Using the Dyson collective field action, 
we show that the free energy obeys the Seiberg-Witten relations. As suggested by Mironov and Morozov, the 
quantized differential form can be obtained from the classical one by the action of a differential operator in the 
hypermultiplet masses and the Coulomb branch modulus. 
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1 Introduction 



Nekrasov and Shatashvili proposed in [1] to quantize the correspondence between the Seiberg-Witten low energy 
effective theory of 4d N = 2 supersymmetric gauge theories [2], and integrable models characterized by a complex 
curve, the Hitchin systems [3]. The Yang- Yang function of the quantum integrable model [4] is associated to a (3- 
deformed prepotential given by the Nekrasov instanton partition function of the gauge theory on the $7-background 
with one vanishing deformation parameter €2 = 0. The second deformation parameter e\ is then identified with 
the Planck constant. When this second parameter is also vanishing, the instanton partition function reduces to the 
Seiberg-Witten prepotential of the gauge theory on M 4 [5, 6]. Based on this correspondence, it has been proposed 
in [7] to compute the quantized prepotential perturbatively in h from Bohr-Sommerfeld integrals appearing in the 
semi-classical treatment of the quantum mechanical system. This prepotential was then successfully compared to 
the instanton partition function at first orders. 

On the other hand, the celebrated Alday-Gaiotto-Tashikawa (AGT) conjecture [8] relate the instanton partition 
function of 4d J\f = 2 SUSY gauge theories with a product of SU{2) gauge groups to a correlator of Liouville theory 
defined on a Riemann surface. 1 This conjecture has been checked in a number of specific limiting cases [11, 12], and 
even proven for a single SU{2) gauge group and Nf < 2 flavors [13]. We focus here on the N = 2 gauge theory 
with a SU(2) gauge group and four scalar matter hypermultiplets in the fundamental representation. This theory is 
known to be scale invariant even at the non-perturbative level, and is sometimes referred to as superconformal QCD. 
The instanton contribution to the partition function is related to the Liouville conformal block, and the full partition 
function, integrated over the moduli space, to the Liouville four points function. Inspired by the CFT methods 
developed for matrix models in [14] (see also [15]), Dijkgraaf and Vafa [16] suggested to realize the Liouville 
conformal block as a /3-deformed matrix model, the /3-ensemble. This proposal was soon generalized to higher 
genus and various gauge groups [17, 18] and thoroughly investigated [19, 20, 21]. 

The proposal exposed in [16] reformulates the AGT conjecture as a relation between the /3-ensemble integrals 
and the Nekrasov partition function. This identity has been proved in [22] for the case (3 = 1 (Hermitian matrix 
model) which corresponds to the relation e\ = —62 between the gauge theory deformation parameters. 2 In the large 
N limit, where N is the size of the matrix, the /3-ensemble reproduces the Seiberg-Witten theory [24, 25], and the 
Seiberg-Witten curve in the Gaiotto form [26] is recovered as the spectral curve of the matrix model. 

The instanton partition function with £2 = is believed to satisfy 7i-deformed Seiberg-Witten relations [1, 7]. 
Following the Dijkgraaf- Vafa formulation of the AGT conjecture, the free energy of the /3-ensemble should also 
satisfies the same relations in the large N, or planar, limit. Until now, this claim has only been verified for /3 = 1 [27]. 
One of the main purpose of this paper is to extend this demonstration using the Dyson collective field representation 
of the planar free energy [28]. When e-deformation parameters are turned on, the differential form of Seiberg-Witten 
theory has to be modified. It is largely believed that the quantized differential form is given by the resolvent of the 
/3-ensemble as H(z)dz [29]. Our study of the collective field action shows that this statement is slightly incorrect, 
and one has to withdraw a regular function from H (z) to recover the differential form di](z)dz. This regular term is 
usually discarded from the cycle integrals of the Seiberg-Witten relations as a total derivative. However, it contains 
poles at the branch points which gives a non-zero contribution, as already noticed in [30]. Such poles do not pertain 
to the singularities of the logarithmic potential but are a general feature of the WKB expansion. The discrepancy 
between the resolvent and the differential form is vanishing at h = 0, it is related to the cut-off term of the Dyson 
action. Moreover, the resolvent H(z) is the classical U{1) current in the free field representation of the /3-ensemble, 
and the relation (2.36) between H(z) and rj{z) is the transformation law of this current under the conformal mapping 
r] — > z(t]). Thus, rj(z) can be interpreted as the change of coordinates which trivializes this classical current. 

Superconformal QCD with SU (2) gauge group is characterized by a four punctured Riemann sphere in the M- 
branes construction [26]. In this simple case, the associated quantum Hitchin system takes the form of a Schrodinger 

'The correspondence was later generalized to SU (n) gauge groups and Toda field theory [9, 10]. 
2 This demonstration has been generalized to an SU(n) gauge group in [23]. 



1 



equation [31, 32]. The quantized differential form can be defined from the solution of this Schrodinger equation, 
and computed explicitly in the WKB expansion. In [7] was suggested the existence of a differential operator in the 
parameters of the Seiberg-Witten curve (mass of the hypermultiplets, Coulomb branch modulus, gauge coupling) 
generating the higher order terms in the ^-expansion of the deformed differential dr](z)dz from the /i-classical 
differential. The existence of such a structure has been observed at first order in h in [31], and will be demonstrated 
to hold at all orders. 

The restriction to the gauge theory with four massive hypermultiplets here is mostly for simplicity, and we 
believe our results extend to theories with a smaller number of flavors. These theories may indeed be obtained by 
sending some of the hypermultiplet masses to infinity, a procedure which has an equivalent in Liouville theory and 
/3-ensemble [33]. We should more particularly emphasize that the demonstration of the planar free energy subjection 
to the Seiberg Witten relations does not depend on the explicit form of the potential. Extending the proof to higher 
orders in the large N, or topological, expansion is a challenging issue that may be addressed using the topological 
recursion tools. This matrix model technique recently extended to the /^-ensembles [34, 35, 36] allows to compute 
the full free energy of the model which is identified with the (ei, e2)-deformed prepotential. Whether this quantity 
satisfies the Seiberg-Witten relations is still an open question. 

It would be interesting to consider further generalizations to higher gauge groups and quiver /3-ensembles. An- 
other promising direction is the relation with the direct formulation of the Nekrasov partition function as a (3- 
ensemble [37]. Finally, one of the most intriguing open questions remains to find a connection with the deformed 
Seiberg-Witten curve introduced in [38], and derived with the method of Young tableaux profile density developed 
in [6] . To do so, we need a definition of the quantized curve which is independent of the quantization scheme. This 
may be given by the underlying integrable models [7, 39, 40], or could also arise from a geometrical point of view 
[41]. 

2 Planar limit of the /3-ensemble 
2.1 Definition of the model 

To formulate the AGT conjecture on the /3-ensemble partition function, we start from the Liouville correlator of four 
vertex operators inserted at the points qj = 0, 1, q, oo and with respective charge aj = ao, «i, 012, cw We further 
assume that the charges satisfy the neutrality condition 

a f + Nb = Q L , (2.1) 

/=0,l,2,oo 

where N is a positive integer and Ql = b + 1/6 is the Liouville background charge. Following [8], we shift 
two of the charges in order to absorb the background charge Ql, ao — > ao + Ql/2 and — > a^ + Ql/2. 
This manipulation corresponds to concentrate the curvature of the sphere equally at zero and infinity. Since we 
are interested in a semi-classical expansion of heavy fields [42], we introduce a scale g and the rescaled quantities 
rrif = —2igctf, h = —igQh? This scale g plays the role of the Planck constant for Liouville theory, and we will 
later consider the expansion g — > 0, N — > 00 at fixed gN, h and m/. In these new variables, the condition (2.1) is 
rewritten _ 

m f + 2y^gN = 0, (with b = iy 7 /?). (2.2) 

/=0,l,2,oo 

After integration over the zero mode, the Liouville four points correlation function exhibits a pole when the rela- 
tion (2.1) is satisfied. The residue can be expressed as a free field correlator with N screening charges inserted [43]. 
The conformal blocks entering in the decomposition of such correlators have a well known integral representation 

3 Our case differ from [12] where the AGT conjecture has been proved in the semi-classical limit c — > 00 (or e 2 = 0) but only for light 
fields, i.e. mf = 0(^2). 
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which was studied by Dotsenko and Fateev [44]. According to the proposal of Dijkgraaf and Vafa [16], further 
studied in [19, 20, 21], these integrals can be written as a /3-ensemble partition function. The crucial issue is the 
specification of the integration contours for the screening operators, it has been resolved in the perspective of the 
AGT conjecture. This choice is believed to be equivalent to the quasiclassical approach to /3-ensembles we use in 
this paper, as shown at first orders in [45]. 

Using the /3-ensemble representation for the Liouville conformal block, the AGT conjecture reads 

%ens(a, rm, 9, P) = q -^-m^ + A K m -K)W (1 _ ^f™,-^)^-^)/^^^ ^ ^ ^ 

up to a (/-independent constant factor that can be found in [20] . 4 On the RHS, ^Nekrasov refers to the full gauge 
theory partition function, including classical, one-loop and instantons contributions. It depends on the UV gauge 
coupling r through the parameter q = e 2t7TT identified with the cross-ratio of the conformal block. On the LHS, 
the /3-ensemble partition function Zp. &D& is a multiple integral over N variables, the eigenvalues Xj, with a potential 
V(X) and a Van der Monde determinant at the power 2/3, 

r N %/? N 

%ens(a, m,9, P) = / 11 e~ y(A/) dA 7 \{ (A/ - Xjf . (2.4) 
J i=i 

iit j 

When /3 = 1, we recognize the eigenvalue integrals appearing in the Hermitian matrix model after integrating out 
the angular degrees of freedom, and the model (2.4) is sometimes also referred as a "/3-deformed" matrix model. 
The "Penner type" potential V contains logarithmic terms taken at the finite insertion points of the vertex operators, 

V(z) = (mo + h) log z + mi log(z — 1) + m 2 log(z — q). (2.5) 

This potential is independent of the parameter m^ which appears in (2.4) only through the matrix size N following 
the relation (2.2). The masses \ii of the four matter hypermultiplets are related to the charge of the vertex operators 
through 

fJ , 1 = m 1 +m 00 , fi 3 = m 2 + m , 

(2.6) 

u, 2 =m 1 -m 00 , m = m 2 - m , 

and with a slight abuse of terminology we call the parameters mj "masses" associated to the insertion points qf = 
0, 1, q, oo. The equivariant deformation parameters e\ and e 2 are connected to the /3-ensemble quantities as 

9 2 = -\^2, /3 = -| = -6 2 , ft= JL(0-l) = I( ei + e 2 ). (2.7) 

Finally, the Coulomb branch vev of the gauge theory a is related to the momentum a of the conformal block in the 
intermediate channel as a — h = 2iga. In the quasiclassical approach, the eigenvalues of the /3-ensemble condense 
on a continuum of saddle points neighboring the minima of the potential, leading to an eigenvalue density supported 
on several intervals. The potential (2.5) defines a model with two such intervals. These are the two branch cuts of the 
resolvent, the Cauchy transform of the eigenvalue density. The filling fraction a describes the ratio of the eigenvalues 
that occupy one of these branch cuts. More precisely, this filling fraction condition will be later interpreted as one of 
the deformed Seiberg-Witten relations. 



4 We consider here the /3-ensemble used in [24, 30]. To obtain the correct powers of q and q — 1 from [8], one should take into account 
the flip a — > —a, mo — > -mo, m x — > —moo, mi — > 2h — mi and mi — > 2h — m,2 that leaves the SU(2) instanton partition function 
invariant. This flip corresponds to the action of the Weyl group of SU(2) gauge and flavor symmetry, and to the reflexion a — > Ql — a of 
Liouville vertex operators. 
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The deformation parameters have dimension of a mass, and the Coulomb branch vev a inherits the dimension of 
the scalar field in the gauge multiplet. The conformal invariance of the gauge theory implies the invariance under the 
rescaling m — > 7/Xj, a — > 7a and ei,2 — > 761,2- 5 This transformation allows to fix one of the e-parameters to one, so 
that the theory depends only on the ratio /3. Similarly, the Liouville theory only depends on the parameter b = iy 7 /? 
and the scale g has been introduced by hand. After a proper rescaling, /3-ensemble correlators only depend on /3 
through the parameter h. 6 Again, the parameters (g, h) are not independent but obey the invariance under rescaling 
g — > 75, h — > jh together with my — > jirif and a — > 7a. 

In the large N limit of the /3-ensemble, we send g — > so that the fixed product gN obeys the relation (2.2), and 
the large N topological expansion is an expansion at small g. There are two possibilities to define this expansion, 
considering either /3 or h fixed. Due to the scale invariance mentioned above, it seems more natural in the context of 
AGT to fix (3 so that h = 0(g). However, it reveals more convenient here to keep the two variables (g, h) independent 
and take the double expansion of the /3-ensemble quantities. In this paper, we should focus on the zeroth order in 
g, but all orders in h, which we refer as the "planar limit". This limit is equivalent to the Nekrasov-Shatashvili 
limit [1] where one of the deformation parameters is set to zero, while the second one is identified with the Planck 
constant 2h. The remaining h expansion of planar quantities is interpreted as a WKB expansion. The planar limit 
also coincides with the semi-classical limit of Liouville theory in which c — > 00. 

Another interesting limit is h — > with fixed g, i.e. (3 = 1, where the /3-ensemble is identified with an Hermitian 
matrix model. It corresponds to set ei = —e 2 for the $7-background, and gives the c — > 1 limit of Liouville theory 
(Ql — > 0). In this limit, the AGT relation (2.3) between the instanton partition function and the /3-ensemble has 
been demonstrated [22]. Combining both limiting cases, we reach the planar limit of an Hermitian matrix model at 
g = h = 0. It corresponds to consider the gauge theory on the flat R 4 background. On the Liouville side of the 
correspondence, we get the semi-classical limit of a c = 1 free boson. 

We define the planar free energy and the deformed Seiberg-Witten prepotential as 

Tp = lim g 2 log Zp, T S w = - lim eie 2 log ^Nekrasov, (2.8) 
and in the planar limit the AGT relation (2.3) reduces to 

J"sw = 47^ + ((mo - m 2 ) 2 + 4k{m - K)) \ogq- 2(mi - 2h)(m 2 - 2h) \og(q - 1) + C (2.9) 

with ei = 2h and a ^-independent constant C. This constant is also believed to be independent of a. The relation 
(2.9) was shown to hold up to the order 0(h, q 4 ) at the special point mo = m^ = 0, mi = m-2 in [30] and this 
result has been extended very recently to the order 0(h 2 ) [46]. 

2.2 Loop equation and Schrodinger equation 

Over the past fifty years, matrix models have been thoroughly investigated and various methods have been developed 
to solve them. Unfortunately, only a part of these techniques have been successfully extended to /3-ensembles [47]. 
Nonetheless, a /3-deformed version of the topological recursion has been developed in [34, 35]. These algebraic 
geometry methods provide a formal recursive solution to the loop equations that still needs to be analyzed in the 
framework of the AGT conjecture. In this context, we believe that the methods and ansatz used here in the planar 
limit have a straightforward generalization to higher orders. 

Loop equations are derived by exploiting the invariance of the measure under a change of variables. The simplest 
equation is obtained from the shift 

A/^A/ + ^— , (2.10) 
z — Xi 

The gauge theory side scale invariance should not be confused with the Liouville theory scale invariance which acts on the insertion 
position of the vertex operators and has already been used to fix these positions at 0, 1, q, 00. 

6 This property can be used to fix (3 upon a change in the matrix size which controls the value of g. At /3 = 1 /2 (or /3 = 2) the partition 
function is identified with the eigenvalue integrals of an orthogonal (or simplectic) matrix model. The limit /3 — ¥ 1 with h fixed corresponds 
to send g to infinity, or the matrix size toward zero. 
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at the order O(e) and has been derived in [25] for the Penner type model (2.5). It writes at the first order in g, 




W(z) 2 + hW'(z) + V'(z)W(z) - U{z) = 0, (2.11) 
where we defined the planar resolvent W and the auxiliary correlator U as 

^w-js^x:^). ^)=iim^E( n 1:r; (Af> )- 

The quantity f7(z) is a sum of simple poles at z = qi where go = 0, q\ = 1, 92 = <Z, with residues proportional 
to W(qi). These three quantities are not independent but are related through two loop equations derived from the 
change of variable Xj — > Xj + e and A/ — > Xi + eXi respectively, 

2 2 

(m + S i0 h)W( qi ) = 0, 4^%(m; + <5 i0 ft)W(%) = m 2 OQ - ( J^m; ) . (2.13) 

i=0 i=0 

These relations can also be obtained from the asymptotic of (2.11) at large z, noting that W(z) ~ yf]3gN/z. They 
allow to express U(z) in terms of the mass parameters, the cross-ratio q, and the quantity u = qni2W(q). An 
interesting feature of the Penner-type potential (2.5) is that the resolvent evaluated at the insertion point z = q is 
simply the derivative of the free energy with respect to this parameter. Combining this property with the planar AGT 
relation (2.9), we may write 

q^^ =Au + {m G -m 2 ) 2 + m{m Q -h)-2{m 1 -2h){rn 2 -2h)-^—, u = q^. (2.14) 
oq q — 1 oq 

The comparison between the matrix model at (3 = 1 and the Seiberg-Witten theory made in [25] relates the quantity 
u, up to a mass-dependent translation, to the Coulomb branch parameter (tr <I> 2 ) where <I> is the scalar field in the 
M = 2 gauge multiplet. 7 In the Seiberg-Witten theory, it is usual to consider u as an independent parameter. On 
the other hand, u being proportional to the planar resolvent taken at z = q, it has a natural expansion in h. The 
relation between these two points of view, u fixed (a function of u and h) and u function of a (fixed) and h, will be 
investigated in the last section, and for now we adopt the Seiberg-Witten perspective. In addition, to simplify the 
quantities, we shift the variable u into 

E = 4u + 2m2(mo + 2h) + 2m\m,2 — — — vn^ — m\. (2.16) 

With a slight abuse of terminology, we will refer to E as the Coulomb modulus parameter of the gauge theory. 
Introducing 

H{z) = W{z) + l -V'{z), T{z) = U(z) + \V(zf + \v\z\ (2.17) 
the planar loop equation (2.11) takes the form of a Riccati equation, 

H{z) 2 + hdH{z) =T(z). (2.18) 



7 The precise identification with the results of [25] involves the change of notations from here to there: mo — > mi, ra\ — > mi, m.2 — > m.3, 
moo — > m and u is related to their quantity U as 

(q + 1)U = —4u(q — 1) + qrrv^ — (2q — l)mj — qm\ + m% — Amimi — 2(q — l)mom2 (2.15) 

in our notations. 
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The meromorphic function T(z) has double poles at z = qi and behaves as 0(1/ z 2 ) at infinity. It has a linear 
dependence in E, 

T(z) = + {q 2 ~ l)E , (2.19) 

where p±(z) is a polynomial of degree four in z which is quadratic in h, it can be expanded as 

P4 (Z) = pf\z) + + H 2 pf\z), 

( O™, 2 i Q,™ 2 , __2\ 2 o™2 | 2 0^,2^ 3 



Ap\ >(z) = rn^z 4 + (g(-2m^ + + mg) - - + m\- 2ml) 



z 



+ (9 2 ("4 - 2ml - mg) + q(2m 2 00 - 2m\ - 2m\ + 2m 2 ) + 2m^ + 2m 2 ) z 2 {2 .20) 
+ (q 2 (-m 2 OQ + 2m 2 , + mf) - q(m\ + 3m 2 ,)) z + m\q 2 



pf(z) = -(J- + £-) = -(* - D 2 (, - q?. 

It is well known that the Riccati equation (2.18) can be cast as a Schrodinger equation for the wave function tp 
defined as H = hd log ip, 

( h 2 d 2 - ) iP(z) = {q 2 ~ l)E Mz). (2.21) 

V y\U(z-<h?) 4n 2 =0 (^-%) 

This one dimensional quantum mechanical problem is the integrable system associated to the Riemann sphere with 
four punctures at z = qj [31, 48]. Up to the pole factor, the parameter E is seen as the energy of the system, whereas 
the ratio in the LHS is identified with the potential. In Liouville theory, the wave function ip(z) is constructed as the 
insertion of the degenerate operator 4>2,i(z) within the four point correlator. The Schrodinger equation (2.21) is a 
consequence of the null state condition of level two obeyed by this operator. In the /3-ensemble, the wave function 
is realized as the planar limit of the following correlator, 



*!>(z) 



e^W lim (j[(z -X^^. (2.22) 



This quantity is associated through the AGT correspondence to the vev of a surface operator in the gauge theory 
[49]. 

The appearance of the Schrodinger equation is at the origin of the "quantum spectral curve" notion defined for 
the /3-ensemble in [36]. In this context, two non-commuting variable x and y satisfying [y,x] = hare introduced. 
The spectral curve is the rational fraction £ (x, y) = y 2 — T(x) and act on the wave functions ip(x) as 

£(x,y)il>(x) = 0, V = ^q- x - (2-23) 

At finite h, the branch cuts of the resolvent disappear, they are replaced by half-line of accumulation of zeros for 
tp. Here, we take a different approach and keep h infinitesimal, using the WKB expansion to solve the Schrodinger 
equation. At each order in h the resolvent keeps the same branch cuts, which is suited for a comparison with a 
/?-deformed Seiberg-Witten theory. 

2.3 Classical (or Hermitian) limit 

The WKB solution consists in expanding the functions in h to solve the problem recursively at each order, 

oo 2 

H{z) = Y i h n H n (z), T(z) = Y J h n Tn(z). (2.24) 

n=0 n=0 
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Figure 1: Definition of the A and B cycles on the spectral curve. 



At the zeroth order, we recover the planar limit of an Hermitian matrix model. The equation for the critical resolvent 
Hq{z) defines the spectral curve (Ho, z) G C x C of the model, 

2 

H {zf= P f Z) p 4{z)=p f) {z) + h q - 1)E f[(z- qi ). (2.25) 

lL=o ( z - aA 1 



i=0 



In [24, 25], this curve has been identified with the Seiberg-Witten curve [2] of the SU(2) superconformal gauge 
theory on M 4 in the Gaiotto form [26]. It describes a torus with singular points at z = qf which is a double cover of 
the four punctured Riemann sphere. In the complex z-plane, the function Hq(z) has four branch points denoted z\ 

tO Zi, 

2 4 

P,(z) = r ^fH(z-z a ). (2.26) 

a=l 

These branch points are functions of the masses "i/e{o,i,2,oo}> the Coulomb branch modulus E, and the gauge 
coupling q. The two dual cycles A and B on the torus will be defined as in the figure 1 . 

In [25], the Seiberg-Witten differential has been identified with dSo = 2Hq(z)cIz, and the filling fraction condi- 
tion provides the first Seiberg-Witten relation [2], 



-(0) 

°< £ > = ti? W {E) = t dS °- 12271 



dSo dF 3W 



The planar free energy of the hermitian matrix model (/3 = 1) was shown to obey the second Seiberg-Witten relation 
(2.27) in [27]. It has been successfully compared to the limit F^ = lim ei ^o -^sw of the Nekrasov partition at the 
point mo = moo = 0, mi = m-2 [24]. However, instead of performing the B-cycle integral, the g-derivative identity 
(2.14) has been used, leaving the possibility to add a ^-independent constant to the prepotential. The classical 
prepotential F^ can be recovered upon integration of the second Seiberg-Witten relation (2.27), after inversion of 
the first relation into E(a). 

2.4 WKB expansion 

Pluging the WKB expansion (2.24) in the Riccati equation (2. 18), we obtain the recursion relation 

n— 1 

2HqHu = T n — dH n -i — ^2 HkH n _k, (2.28) 

k=l 

where T n vanishes for n > 2. This equation has been solved up to n = 2, and used to recover the free energy 
at the order 0(h 2 ) [30, 31, 46]. It is possible to go further by using the ansatz H n = Hn^ + Hn^ Hq with two 

(r) (s) 

meromorphic functions Hn and H n . The regular term of this decomposition arises naturally from the Riccati 
equation and can be observed at the first orders. The second term assumes that the function has the same branch 
points as the ^-classical current Hq(z). This is in the spirit of [7] where the quantum prepotential is obtain from the 
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cycle integral of a deformed differential form, the underlying curve remaining classical. The opportunity to keep the 
classical spectral curve is due to the possibility of expanding in h the branch points z a (h) of a /3-deformed curve, 
\J z — z a (K) being at each order 0(h n ) a meromorphic function with a pole of degree n at z = z a {0) times the 
"classical" square root \J z — z a (0). In this setting, branch points of negative order arise naturally. We should stress 

(r) (s) 

that these poles of Hn and H n are not related to the singularities of the Penner type potential (2.5) but are a 
general property of the WKB expansion. 

The recursion relation (2.28) splits into two relations for regular and singular parts respectively, 

n-1 



k=2 

rp -, / / / ' : \ , II — - 

oo-(s) _ J « STtt(s)tt{s) 1 ft ( I 1 rijHffH 

2^ -3b-Xj ff * ^ 



k ~ 2 , , (2-29) 

n- 1 -. / rr(r) \ , n-2 



where n > 2 and we used the expression of h\ = — (1/2)8Hq/ Hq. 

(r) (s) 

We now analyze the singularities of the functions H n and H n , we expect poles located at the branch points 
z = z a and at the singularities z = qi of the spectral curve. We first notice that the ratios T n /To have only simple 
poles at z = z a , and no poles at z = qu they behave at most as a constant when z — > oo. Next, we look at the first 
terms of the recursion (2.29), 

1 y> 2j^z-q % A^z-Za' iy) 2T (z) 

(r) 

These two functions have simple poles at the branch points z = z a , and H\ has additional poles at z = q{. At 
infinity, behaves as l/2z and as 0(1/ z). 

(r) (s) 

From the recursion relation (2.29), it is possible to show that the functions Hn and H n with n > 2 have no 
poles at z = qi. We further see that H^ is at most 0(1/ z 2 ) at infinity, and H n s ^ at most constant. These constants 
are actually vanishing due to the constraint W(z) ~ yf]3gN/z at z — > oo which imposes H n (z) ~ 0(1/ z" 2 ) for 
n > 2, and consequently H^^z) ~ 0(1/ z). Finally, we denote a$ and a$ the order of the pole at z = z a for 
the functions H^ and H n s \ The recursion (2.29) suggests affl = a^li + 1 an d a[ s ' = a^]_ x + 2 which, together 

with the initial value af^ = df* = 1, is solved by = |_(3n — l)/2j and q^- 1 = [3n/2j where [a?J denotes the 
lower integer bound of x. Coming back to the two equations (2.29), we check that indeed the orders of the poles at 

(r) (s) 

z = z a in each terms of the LHS have upper bound or n and a n respectively. These analytical properties allow to 

(r) (s) 

write the functions H n and Hn as 

p=l ce=l p=l a=l 

with moments /x and /i independent of These expressions are convenient to compute the A and B cycles integrals 
of H(z) which can now be formulated as hypergeometric functions. They also provide a good starting point for 
the star operator construction in the CFT methods developed in [15] since they exhibit the singular expansion of the 
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classical current H(z) at the branch points. The moments fj, and jl obey a recursion relation inherited from (2.29), 



(r) (s) (r) (s) 

-1 /<*fc ' a n-k <*k a n-k 



,(«) = P ~ 1 ~(n-l) 

fc=2 \ <J=1 r=l <j=p r=l 



_ r (n~i) _v VV 4- V V Vf-iv-*?^- 1 z p-i- r u (k) r X n - k ) 

Mp,a — 2 Z^ I Z^ Z_^ ^g,a^r,a p _ g _ r -f ^ ^ ^ \ L ) ^q+r-p- l Z af3 ^q^rfi 



and, 



+ EE E(- 1 ) p " r ^ 1 - P -i^7" r ^>fe- fe) 

<j=l r=p p^ct / 



(r) M 

(n) ^ ^ n-1 4_ Vfe _ 

,W ,,(n-fc) -q 



(2.32) 



m co I l^a z aP k=lq =l r =l 0^ 



'r,7 ' Z «/3 Z /?7 

n-1 /"fc Q „-fc «L <-fc \ 

-E EE /42/4ir fe) + 2 y: E E(- i r^S-P-i^ 9 " r «^7 fe) 

fc=l V q=l r=l q=p r=l P^a / 

(r) (r) (r) (r) ( r ) (r) 

a „-l "n-1 n-l a fc + a n-k / 4 Q n-fc 

« Z^ Z^ « ^g,/? Z "/3 Z^ Z_-/ I Z^ Z^ ^q,a^r,a Op+m-q-r-1 

q=p—l q=l f3^a k=l m=0 \ q=l r=l 

(r) (r) 

a; ; a . \ 
o i\p+m-q r r-l p+m-q-r-1 (fc) (ra-fc) \ 

Z Z^ 2^ 2^ ' ^g+r-p-m^a/? N^r,^ I 

g=p+m— 1 r=l /3^a / 

(2.33) 

where we denoted z a( g = z a — zp, the binomial coefficient Cp, and 

oo 

:. ^oW- 2 = E( z -^) n " lr « n) - ^ 234 ) 



r 



(n) = £ (f ~ " d * 

F U) 2 2ivr' 



2.5 Formal resummation 

The previous /i-expansion for singular and regular parts of H can be formally resummed into the functions H ^ and 
H( s \ These two functions are related to each other through 

ffM(z) = -(h/2)d\og(H^(z)H (z)), (2.35) 

which is a resummed version of the first recursion relation in (2.29). We further introduce the function rj(z) defined 
up to an additive constant by H^Hq = drj, so that = — {h / 2) d log drj. Reassembling singular and regular 
parts, we write the relation between rf and H as 

H(z) = dr,{z) - ^d\ogdr,{z). (2.36) 
Inserting this definition into the Riccati equation, we find an equation over rj which takes the form 

T(z) = (d V (z)) 2 -^{ V (z),z}, (2.37) 
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where {rj, z} denotes the Schwartz derivative. This expression arises naturally in the Schrodinger equation context, 8 
and appears in the accessory parameters problem (see for instance [50]). Indeed, combining the relation between h, 
H and ip, and given that rj and — rj satisfy the same equation (2.37), the two independent solutions of the Schrodinger 
equation write 

/i" p ±/?, _1 »;(z) 
2^g> (2-38) 

where the constant factors are fixed by taking the Wronskian ip~dip + — tfj+dip" equal to one. At the order 0(h), 
dr] is equal to the classical momentum. The ratio of the two solutions, r = t/j~ /ip + is known to satisfy T = 
— (h 2 /2){r, z}, replacing r = e~~ 2h v we recover the equation (2.37). 

There is another interesting interpretation for the expressions (2.36) and (2.37). In the CFT formulation of matrix 
models [14, 15] extended to /3-ensembles [16, 17], the partition function Zp reproduces the correlator of a free field 
with a background charge proportional to h. Up to a Wick rotation, this field is the holomorphic part of the Liouville 
field after integration over the zero mode. In this framework, the planar limit is the classical limit of this Coulomb gas 
field, and H(z) is interpreted as the ^-classical U(l) current H = d(f>. Similarly, T(z) corresponds to the (/-classical 
stress-energy tensor, and the Riccati equation (2.18) simply express the standard definition T = (dcf)) 2 + hd 2 cf). 
The relation (2.37) between T(z) and rj(z) can be seen as the transformation of a constant stress-energy tensor 
under the inverse mapping rj — > z(r]). 9 When the background charge is non-zero, the U (1) current is no longer a 
primary field, and the relation (2.36) is the transformation of a constant current under the same conformal mapping 
rj — > z(rj). Thus, in the coordinate rj, both classical stress-energy tensor and U (1) current are equal to one. Solving 
the Schrodinger equation in this coordinate, we get 

h2 ^(v) = Hv) ^(V) = yfe ±fi "\ (2-39) 

which indeed gives H(r)) = hd v log ^ + (??) = 1. The transformation law of the wave functions V ,=t can be deduced 
from their relation to the classical current, they transform as primary fields of dimension —1/2 under rj — > z(rf), 

Y^r?) -> (driiz))- 1 / 2 ^^)), (2.40) 

and we recover the expression (2.38) from (2.39). 



3 Deformed Seiberg-Witten theory 

3.1 Planar free energy and Seiberg-Witten relations 

We now show that the planar free energy Fp of the /3-ensemble obeys the Seiberg-Witten relations with a deformed 
differential form dS = 2di]dz given by only the singular part of the resolvent, 




(3.1) 



with the cycles A and B previously defined on the classical spectral curve. Let us emphasize that we consider here 
the formal resumation in h, and these relations, together with all the identities below, only hold order by order in the 
h expansion. In the limit h — > 0, drj reduces to Ho and we recover the standard Seiberg-Witten relations (2.27). 
To analyze the first relation in (3.1), we decompose H into regular and singular parts as in (2.36), 

/ H(z)dz = I dr,dz - - I 9 ^°,, dz - - I d\og(H^(z))dz. (3.2) 
J A J A 2 J A Hq(z) 2 J A 

8 We thank V. Pasquier for suggesting this connection. 

9 In the context of Hitchin systems, this is the transformation law of a projective connection [39]. 
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The last term in the RHS is proportional to the monodromy of \og(H^ (z)) around the cycle A. Since the expansion 
of the logarithm in h produces only regular terms at each order, it is vanishing. The second term in the RHS of (3.2) 
exactly compensates the shift of — h for the filling fraction condition observed in [30], 

I 2H(z)— =a-h. (3.3) 
J A 2m 

This relation is believed to be exact at the order 0(h 2 ) [46]. 

To demonstrate the second relation (3.1), we employ the collective field formulation of the /3-ensemble at large 
N derived by Dyson [28]. A similar calculation has been done in [27] for h = and we extend it here to the f3- 
deformation. Closely related considerations are found in the study of complex /3-ensembles in which the eigenvalues 
are integrated over the full complex plane, leading to a bidimensional Dyson gas [51]. 

We first review the derivation of the Dyson action [28]. In the large N limit, the /3-ensemble partition function 
can be written as a path integral over the density of eigenvalues p{x) seen as a collective field, 

Z P = J D[p]e^ A[p \ p{x) = g^l^y{x-\ I )\. (3.4) 

Since we work perturbatively in h, the density at each order has a finite support F corresponding to the branch cuts 
of the classical resolvent Hq{z). Again, this is no longer true at finite h where the eigenvalues remain isolated. The 
action appearing in (3.4) is obtained from the effective action for the eigenvalues A/ deduced from the definition 
(2.4), 

i m N N N 

^4eff[A,] = ^Vv(A/) + /3 V log|A J -A J |-^VlogZ(A / ), (3.5) 
y y 1=1 I,J=l 1=1 

where we introduced a short distance cut-off l(x) to regularize the sum of logarithms at coincident indices I = J. 
This cut-off is interpreted as the mean distance between the eigenvalues at the point x, it also arises in the Jacobian 
of the transformation, 

N N 

H d\j = N\J[p]D[p], J[p] oc Y[ /(A/). (3.6) 

7=1 1=1 

Combining these two contributions, and introducing the density within the effective action, we derive the following 
collective field action 

A[p] = / V(, W *>fa + / log |, - yWiMy^y - h J „(*) log «(«)* + est. (3.7, 

Since a "single" eigenvalue roughly occupies a segment of size proportional to the inverse of the local density, 
l(x) oc l/p(x) and the classical action reads 

A[p] = jv(x)p(x)dx + j log \x - y\p(x)p(y)dxdy + hj p(x) logp(x)dx, (3.8) 

up to a constant set to zero by renormalization of Zp. 

To confirm the expression (3.8) of the collective field action, we should analyze the equation of motion, 



oxop(x) J-p x — y pyx 



where the integral is considered as a principal value. The density solving this equation of motion is the (/-classical, 
or planar, density, i.e. the limit of the density p{x) defined in (3.4) as g — > 0. Since we work only at the planar level, 
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and the full density p(x) will play no role below, we will keep the notation p(x) now referring to the planar density. 
It is related to the planar resolvent (2.12) by a Cauchy transform 



W(z) = [ W(x±iO) = T^v^p(x) + / x G F. (3.10) 

Jr z ~ V Jr x ~ y 

Taking the discontinuity of the loop equation (2.1 1) over the branch cut T, we recover the equation of motion (3.9). 
In particular, the derivative of the resolvent which characterizes the /3-deformation is at the origin of the logarithmic 
derivative term d log p in (3.9). 

The classical spectral curve (2.25) describes a two cuts solution: the density support T is made of two disjoint 
intervals T = T_ U T + with T_ = [z\, Z2] and r + = [zs, Z4]. The moduli of the solution are fixed by imposing the 
two filling fraction conditions 

jf p(x)dx=^(a-H), p (x)dx = y/pgN-±(a-H). (3.11) 

The first identity is equivalent to the Seiberg-Witten relation (3.1) involving the „4-cycle integral which circles the 
branch cut T_. The second relation reflects the normalization of the density to yffigN. 

The planar free energy is the collective field action evaluated for the planar density that solves the equation of 
motion (3.9). Taking the a-derivative of the expression (3.8), we get 

l^ = /^r c(a;) ' withC ( x ) = v ( x ) + 2 j r lo s\ x -y\p(y)dy + nio g p(x). (3.12) 

The equation of motion (3.9) implies that C{x) is a constant on each support T± that we denote C±. In the standard 
approach to matrix models, these two constants are supposed to be equal,which provides an extra requirement on 
p [52]. As mentioned in [27], this is no longer the case in the Seiberg-Witten theory. The derivative of the filling 
fraction conditions (3.11) with respect to a supplies two expressions that can be used to rewrite the free energy 
derivative as the difference of the two constants. 

°g-\tf!-C + ). (3.13) 

To rewrite this equation as a cycle integral, we need to extend the function C(x) to the whole complex plane, i.e. we 
are looking for a function C(z) with a branch cut on V where it satisfies C(x + iO) + C(x — iO) = 2C(x). This 
function can be chosen up to a constant term since we are only interested in the difference C_ — C+. Once it has 
been properly defined, it is possible to write 

dTp 1 f Z3 dC 



i r Z3 dc i r 

da =--^- d z=-J g SC(^, (3.14) 

provided that dC has no regular part. When (5 = 1, the function C(z) is the classical bosonic field appearing in the 
CFT methods [14, 15, 16, 17], C(z) = 2<p{z) with H(z) = d<j>(z). When h ^ 0, we need an additional contribution 
C(z) = 2(f)(z) + hL(z) because of the logarithmic term. This function L(z) have a branch cut on T where it verifies 

L(x + iO) + L(x - iO) = 21ogp(x). (3.15) 

To identify the function L(z), we notice that since drj is by definition the singular part of H, it satisfies dr](x ± iO) = 
=Fi7r/?(x) for x G T. It implies 

log(dr)(x + iO)) + log(dr](x - iO)) = log(^ 2 p(^) 2 ), (3-16) 
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which leads to identify L with \ogdri up to an irrelevant additive constant. We conclude that C(z) = 2rj(z), and 
(3.14) provides the second Seiberg-Witten relation (3.1). The difference = (H/2)dL between the two functions 
H{z) and drj(z) is specific to the /3-deformation. It is related to the cut-off part of the Dyson action, i.e. the term 
p log p in (3.8). The ,6-cycle integral of the relation (3. 1) is of course invariant under a regular shift of the differential, 
but again has poles at the branch points which require a special treatment. 

To summarize, we have shown that the Dyson representation of the planar free energy, 

Tfi = jf ^ \V{z) + Mogdr)(z)} d V (z) + jf ^ jf ^ log \z - zf\drt{z)dr,(zf), (3.17) 

provides a new solution to the Seiberg-Witten relations (3.1). The potential V(z) is chosen such that the function 
T(z) defined in (2.17) satisfies the equation (2.37) in which the RHS is entirely defined by the differential form 
dS = 2dr\dz. In (3.17), the contours circle both branch cuts T_ and T+ but exclude the logarithmic branch points 
of the potential. In addition, the double integral in the RHS must be properly regularized at coincident point z = z' . 



3.2 Differential operator 

In their study of the deformed Seiberg-Witten relations for the pure SU{2) M = 2 gauge theory, Mironov and 
Morozov [7] introduced a differential operator in the Coulomb branch modulus and the gauge coupling to derive the 
0{h 2 ) corrections to the quantized differential form dS. A similar operator with derivatives in the parameters of the 
potential was observed in models with Gaussian and cubic potentials [53]. In SU(2) superconformal QCD, such 
a differential operator also appeared at first order in 0(h) [31], with a differentiation with respect to the masses. 
Here we demonstrate the existence of such operators at all order in h. These relations ensures that the quantities 
considered in the last subsection are finite. 

(s) 

The summation of fractions in the expression (2.31) of Hn can be reorganized as 



» 



P =i 



cipZ 3 + b p z 2 + c p z + d p 
Pa{z)p ' 



(3.18) 



A recursion relation on the numerators can be derived directly from (2.29), and would involve Euclidean divisions of 
polynomials. Here we do not need the explicit expression of the coefficients a p , b p , c p and d p which are functions of 
the masses, the Coulomb branch modulus, and the gauge coupling. The degree four polynomial P&(z) decomposes 
as a sum of monomials, 



P 4 (z) 



^z' + hz' + hz' + hz + fo, 



(3.19) 



where the coefficients fj depends on the four masses mj, as well as the parameters q and E. Taking the derivative 
of the spectral curve (2.25) with respect to fj, we easily show that 



dfjWo) Ho{z) - { - l) 2 r[i/2] p,(zy Ho{z) - 



(3.20) 



Using this identity, we rewrite the n-th term in the ^-expansion of the differential drj as a differential operator in fj 
acting on Hq, 



H^Ho 



p-1 



2r[i/2] / d . 



d 



d d d 

W2 +Cp—+d p — 



H 



(3.21) 
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It now remains to change variable from fj to the three masses mo, mi, m<i and the Coulomb branch modulus E. The 
Jacobian matrix of this change of variables is given by 



/ _d_ \ 
/ ah \ 



9 _t 



V 9/o / 



2q 




2 


2 


mi 


(9-1) 2 
2 


«>2 (g-i) 2 

2 


mi 


(?"1) 2 
2 


m 2 (q~l) 2 q 


mi 




2 




2 







2 



(9-1)* 
4(3g-l) 



\m 2 (g-l) 2 g 2 m 1 (q-l) 2 m q 2 



(9-1)* 
8(2g-l) 

(9-l) 3 9 
4(g 3 -3g 2 +8g-4) 

(«-l)V 



/ -2- \ 

9 
dmi 

d 
din 

Km/ 



(3.22) 



Its determinant is non- vanishing as long as the four punctures qj of the Riemann sphere remain separated, i.e. q ^ 
and q / 1. At first order in ft, we recover the result found in [31], 



H[ S) H 



d 



+ 



o 



dm\ dm2 



H . 



(3.23) 



This differential operator does not involve the derivative with respect to the Coulomb parameter E. This simpli- 
fication is due to the previous change of variable (2.16) from u to E. We would have obtained exactly the same 
property if we had chosen instead of E the variable qdJ^sw/dq related to u through (2.14). However, the intermedi- 
ate formulas such as the expression of the polynomials (2.20), or the Jacobian matrix (3.22), would have been more 
complicated. It is not clear priori whether the higher order operators would involve a differentiation with respect to 
E, or if it can again be eliminated after a redefinition of the Coulomb branch parameter. 



3.3 Interpretation of the differential form in Seiberg-Witten theory 

We have shown in subsection 3.1 that the planar free energy obeys the Seiberg-Witten relation with the differential 
dS = 1dr\dz. It remains to see whether this differential is relevant to the Seiberg-Witten theory. In particular, since 
the Seiberg-Witten differential is commonly identified with the full resolvent 2H(z)dz, we should compare here 
H with dm We have seen previously that the „4-cycle integral reproduces the correct result, including the shift of 
— h observed in [30]. The Bohr-Sommerfeld quantization rule applied to the underlying quantum integrable model 
provides a simple argument in favor of replacing H by dm the classical momentum appearing within the cycle 
integrals is equal to drf(z) up to terms of order 0(h 2 ), and to H(z) only up to 0(H). In this context, the shift of —h 
has indeed to be taken into account. 

We now focus on the ,6-cycle integral. As a consequence of the relation (3.23), the planar free energy obeys the 
following identity, 

d ^(1) ( 9 d \ d ro) 



8a P 



8a P 



(3.24) 



dm i dm2 y 

where the superscript pertains to the ft-expansion of Fp. The equality (2.9) between the prepotential and the free 
energy has been checked at first orders in q up to a ^-independent term in [30] . A priori, such a g-independent but 

(r) 

a-dependent term could be related to a non- zero B-cycle integral of H\ , and would break the relation (3.24) for the 
Seiberg-Witten prepotential. In order to discard this term, we have to verify that the prepotential Fsw at the order 
0(q°), i.e. the one-loop contribution, satisfies (3.24). 10 The expression of the one-loop correction to the prepotential 
can be found in [8], 

^SW.l-loop = -^vectorM + ^ -^hyperOi , ±a), -FvectorCa) = 72(2d) + 7 2 (2a + 6l), 

± (3.25) 
^h yP er(^, a) = -72 (a - Hi + ei) - 7 2 (a - H2 + £2) - 72 (-a + M3) - 72 (-a + M4), 

l0 This is trivial for the classical part Jsw.d = a 2 log q. 
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where we separated the N = 2 vector and hypermultiplet contributions, and denoted 72 the limit 62 
logarithm of Barnes' double gamma function r2(x|ei, £2), 

72 (x) = lim eie 2 logr 2 (x|ei,e2). 
The ei -expansion of the function 72 satisfies the remarkable property 

-2 j2 



72 (x) 



of the 



(3.26) 



(3.27) 



Replacing the derivative "/(x) by a differentiation with respect to the m, and then changing variable from /jj to m/, 
it is possible to show that 



hyper 



9 



+ 



3 



<9toi 9m2 



a fixed 



-^hyper(/ i *' a )' 



(3.28) 



In this calculation we considered a fixed, i.e. independent of the variables \ii (or m/) and h. Since at the order 0(e\) 
the contribution of the vector multiplet is vanishing, the previous relation extends to the whole one-loop contribution, 



(i) 



SW,l-loop 



+ 



dmi dm2 



T 



(o) 



a fixed 



SW,l-loop- 



(3.29) 



We have not shown yet that the deformed Seiberg-Witten prepotential at the order O(q ) obeys the relation (3.24). 
There is indeed an important subtlety since we worked with a fixed, and the relation (3.24) holds only for E fixed. It 
is actually possible to rewrite this relation at a fixed, taking care of the reshuffling of the terms in the ^-expansion. 11 
Then, this result dismisses the regular contribution H[ r ^ and confirms the identification of dS = 2dr/dz with the 
Seiberg-Witten differential, up to terms of order 0(h 2 ). 

Another important feature of the Seiberg-Witten differential is the residues at the punctures. The classical dif- 
ferential dSo = 2Ho(z)dz is known to have residues m/ at the z = qf € {0, 1, q, 00}. Both 2H(z)dz and 2r](z)dz 
get ^-corrections that take the form of constant shifts: tuq — > uiq + h and — > m,^ — h for the first one, and 
mi 2 — > fn\,i — h for the later. Up to now, there is no interpretation for such a shift, but it is possible to correct the 
differentials by adding a pole contribution at z = qj that does not modify the A and B cycle integrals. 



4 Discussion: a direct approach to the prepotential 

In the Seiberg-Witten context, E is seen as a parameter independent of h, and a is a function of E and h. This 
is rather unusual from the matrix model perspective where E is a function of a, which is fixed, and h. But both 
approaches are related: E is an intermediate parameter of the Seiberg-Witten theory, it can be taken as any function 
-E(a,fixed> h). When this function is the inverse function of a(E, h) determined by the first Seiberg-Witten relation 
(3.1), we recover the matrix model point of view. 

It is instructive to take this later standpoint, and compose the WKB-expansion of E and H, thus defining a new 
expansion of the classical current [31], 

00 00 
E = J2^En, H(z\E) = J2^H n (z\E ,--- ,E n ), (4.1) 

ri=0 n=0 

"Let a and g be two functions of E, h and some parameter m. Changing variable from E to a, the fi-expansion g(E, ft) = g ( -°\E) + 
hg w (E) + 0(ft 2 ) is reshuffled into g(E(a, ft), ft) = 3 (0) (a) + ftff (1) (a) + 0(ft 2 ). The relation g (1) (E) = d m \ E g {0 \E) is equivalent to 
g^(a) = d m \ a g < ' \a) provided that the second term in the ft-expansion of a(E, ft) obeys = 9 m a (0 \ since g(°\E) = <?' ' (a^ 0) (E)) 
and d m \ a = d m \E - d m a^ {E)\ E d a(a ) ■ 
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where we highlighted the dependence in E. The WKB method presented in the section 2.4 still works in the same 
way, replacing H n by H n , E by Eq and shifting 

T n ^f n = T n + (g - 1 )^ for n > ! (4-2 ) 

4E[Lo (*-*) 

the quantities T n with n > 2 being now non-vanishing. We notice that, i^o(z|i?) = £T(z|-Eo), an d H n still de- 
composes into regular and singular parts which obey similar analytical properties. In particular, we may adapt the 
decomposition (2.31) to H n and write 

^ co = E E - *■)-". ^\ Z ) = - 8 { -^^ n (*-*) + E E - *.)-*■ (4-3) 

p=l a=l i=0 p=l a=l 

Because of the shift in H n , the coefficients ^ip )C i and /x^a depend only on Eq, ■ ■ ■ , E n -\, they obey a slightly 
different recursion relation than (2.32) and (2.33) which is derived in a similar way. The Seiberg-Witten relations 
should also be expanded in h, the zeroth order determines the relation between a and Eq as 

j>^2H (z\E )—=a, (4.4) 

where now a is a fixed parameter. The next orders give the dependence of E u (Eq, ■ ■ ■ , E n -i) which is, up to a q 
integration, the planar free energy of the /^-ensemble, as seen from the relation (2. 14), 12 

p=l a=l v ' 

This method produces directly the prepotential Tp, up to terms independent of q. And again, a differential operator 
may be used to obtain the contour integral from (4.4). 

The previous calculation generates the /3-deformed electric prepotential Fp, but it may also be interesting to 
determine directly the dual magnetic prepotential, related to Fp by a Legendre transform [54]. Replacing the A- 
cycle by £>-cycle integrals in the formulas (4.4) and (4.6), we obtain Eq{cid) and E n (Eo, • • • , E n -±). Finally, we 
may use the identity (2.14) with Tp replaced by and the (/-derivative at fixed a by a derivative at fixed op to 
recover the dual prepotential upon integration. Such a formulation may help to understand the S-duality from the 
/3-ensemble point of view. 
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